Structural change of vortex patterns in anisotropic Bose-Einstein condensates 
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We study the changes in the spatial distribution of vortices in a rotating Bose-Einstein condensate 
due to an increasing eccentricity of the trapping potential. By breaking the rotational symmetry, 
the vortex system undergoes a rich variety of structural changes, including the formation of zig- 
zag and linear configurations. These spatial re- arrangements are well signaled by the change in the 
behavior of the vortex-pattern eigenmodes against the eccentricity parameter. This behaviour allows 
to actively control the distribution of vorticity in many-body systems and opens the possibility to 
study interactions between quantum vortices over a large range of parameters. 
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The superfluid nature of atomic Bose Einstein conden- 
sates (BECs) is one of the most striking manifestations 
of quantum mechanics on a macroscopic scale. Its tell- 
tale sign, the formation of quantized vortices, has been 
extensively studied in recent years and led to significant 
progress in understanding the phenomenon [l|. However, 
to create these topological defects one usually needs to 
break the rotational symmetry of the condensate. This 
requires a significant external disturbance through, for 
example, optical phase imprinting techniques [2[ or stir- 
ring laser fields that allow to excite quadrupole-mode res- 
onances [3j- The latter one is similar to the classical ro- 
tating bucket method that creates a vortex by rotating 
a bucket full of water and was first used in experimen- 
tal studies of superfluid 4 He [J]. The validity of such an 
analogy is certainly limited, given that the non super- 
fluid component in superfluid 4 He experiences friction 
due to its relative motion with respect to the walls. As 
in the classical case, this implies a transfer of energy and 
angular momentum from the walls of the bucket to the 
superfluid part through the non-superfluid component. 
In a trapped BEC the latter is always negligible and the 
transfer has to be made through the excitation of normal 
modes, usually quadrupole ones. A BEC thus reacts to 
a large amount of angular momentum by creating many 
vortices with winding number equal to one [5[, which 
arrange themselves in geometrically defined spatial pat- 
terns. At large vortex density, these structures mimic 
the celebrated Abrikosov lattice [6j. For harmonically 
trapped alkali condensates this was first observed in the 
seminal experiment by Abo-Shaeer et al. where more 
than 100 vortices formed a triangular-shaped lattice with 
a few seconds lifetime. 

Recently, numerical evidence has been provided that 
the vortex pattern of a 2D BEC in an in-plane anisotropic 
rotating trap can undergo structural changes as a func- 
tion of the eccentricity. Specifically, in Ref. [7] it has 
been shown that, for modest changes in the eccentricity, 
an off-line configuration (typical for an Abrikosov lattice) 
can change into a linear one. While this bears analogies 
with the case of ionic crystals [8] , the characterization of 
structural changes in anisotropic and rotating BECs re- 



mains largely unexplored. Most of the existing literature 
focuses on the limit of large numbers of vortices for either 
a symmetric trap [jj Ho) or very high angular frequencies, 
which leads to stripe-shaped vortex patterns jllHUj . Al- 
though the case of medium vorticity has been addressed, 
the role of external forces on the dynamics of the vor- 
tex structures still awaits a systematic approach [l5l-[l8| . 
Yet, understanding how vortices behave under external 
perturbations is a pre-requisite for harnessing the quan- 
tum properties of vortex patterns. Here we present a sig- 
nificant contribution to advance these aims by studying 
the behavior of finite-sized vortex patterns in 2D BECs 
confined within a rotating anisotropic trap. In particular 
we investigate in detail the effects of the eccentricity on 
the spatial distribution of the vortices. By minimizing 
the eccentricity-dependent interaction potential between 
vortices, we show that the vortex configuration 

undergoes structural changes as the eccentricity pa- 
rameter is varied. A hydrodynamical approach to the 
description of the superfluid motion allows us to identify 
the eigenmodes of the vortex-patterns and connect the 
appearance of discontinuities with the transition points 
between different structures. In fact, the modes suggest 
that the change in the equilibrium positions of the vor- 
tices is due to the re- arrangement of the superfluid ve- 
locity field. 



I. VORTEX PATTERN 

We consider the pattern of vortices in the ground state 
of a BEC held in a rotating trapping potential. The 
ground state is found by minimizing the energy func- 
tional [H| 
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where ^ is the normalized order parameter of the conden- 
sate (its dependence on r is omitted for ease of notation), 
V(r) is the trapping potential, m is the atomic mass, N 
is the number of atoms, g = Aith 2 a/m is the inter-atomic 
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interaction energy volume determined by the s-wave scat- 
tering length a, ft is the rotation frequency vector of 
the condensate and L is the angular momentum oper- 
ator. The function \I/ minimizing £ has been studied 
both numerically and analytically under different work- 
ing assumptions such as the Thomas- Fermi (TF) approx- 
imation [19j, the lowest-Landau- level (LLL) approxima- 
tion [20| or the limit of very weak interactions [16]. The 
first usually corresponds to the requirement of a very 
large number of particles, so that the kinetic energy as- 
sociated with VI^atsI (with ^ns representing the non- 
singular part of the order parameter) can be neglected 
in favor of the boson-boson interaction. In the LLL ap- 
proximation, on the other hand, the main contribution to 
the energy stems from the centrifugal term and ^ is well 
described by means of single-particle wave- functions. Fi- 
nally, in the limit of weak interactions the healing length 
becomes large and even under strong rotation only a 
small number of large vortices nucleate [l6| . 

Here we consider a BEC in a harmonic trap rotating 
about its z-axis, which is also the direction of tight- 
confinement, so that ^ can be factorized into an axial 
part (the ground state of a harmonic potential) and an 
in-plane one, ijj(x,y). We call ujj (j=x,y) the trapping 
frequency along axis j of the trap and introduce the ec- 
centricity parameter A = uo y /uo x . 

We are now in a position to minimize £ in the TF limit. 
For a set value of 0<A<1, we call Qn v (X) the minimum 
angular frequency of the trap which allows for N v vor- 
tices in the state which minimizes £?[\I/, \I/*], while is 
the position of the i th vortex in the frame rotating with 



the condensate. By introducing |r, 
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ergy of the vortex pattern can be written as U =Ut 
with [H| 
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Here, po(X)=y/2X/ 7r is the density of the condensate 
at the center of the trapping potential and minimizing 
these energies will determines the positions of vortices. 
In doing this, we will assume that the variations of A 
are accompanied by an adiabatic change of the angu- 
lar frequency so that Qn v (\) < ft <C £1a^+i(A), which 
ensures that the wave-function minimizing the energy 
functional carries exactly N v vortices. The absence 
of a ^-dependent term from the expression of Ut can 
be understood by considering that, in the TF limit, the 
centrifugal force is proportional to the restoring term for 
Q w ^at^(A) (see Eq. (3.4) in Ref. [14]), so that they sum 
up to a quantity which is independent of Q. In order to 
quantitatively asses the deviations of the vortex pattern 
from the Abrikosov-like lattice [7] , we first show how the 
distances of the vortices from the tight trapping direction 
vary against the eccentricity A. Two representative cases 
(N v = 7, 8) of the general dynamics are shown in Fig. [TJ 
the pattern of vortices corresponding to values of A larger 
than a critical threshold Xl (in general a function of N v ) 




FIG. 1. (Color online) Distance Ax of the vortices from the 
soft trapping axis (in units of ^2NgQ J /hjj x ) against the ec- 
centricity of the trap, A. We show the cases of N v = 7 and 8 
(panel (a) and (b) respectively) and plot only the changes in 
positions of four vortices in the lattice (the association with 
the curves is irrelevant) , with the remaining showing analo- 
gous behavior. At A = Xl the vortices suddenly align along 
the y-axis (Ax = 0). We have used a BEC of 10 6 87 Rb atoms 
with a scattering length a — 5.23xl0 -9 mina trap of frequen- 
cies oj z /2tt = 100Hz and y/uj x uj y /2tv = 50Hz (independent of 
A). 
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FIG. 2. (Color online) Phase distribution of the condensate 
carrying a vortex-lattice in the x — y plane with N v = 8 for 
different values of A. The black d ots mark the positions of the 
vortices (in units of ^2NgQ/hoj x ) The appearance of different 
structural vortex patterns is clearly visible. 



abruptly collapses to an all-aligned configuration. 

However, looking at the distance of the vortices from 
the soft axis only gives limited information about the 
actual vortex pattern and we show in Figs. [2] and [3] the 
full position distribution for different numbers of vortices 
and different values of 1>A>Al. Two more structurally 
distinct configurations become evident from this and let 
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FIG. 3. (Color online) Phase distribution of the condensate 
carrying a vortex-lattice in the x — y plane with N v — 7 for 
different values of A. All other values are as in Fig. [2] 



us first consider the case of an even number of vortices 
(shown in Fig. [2]): starting from an Abrikosov-like pat- 
tern at zero eccentricity (A = 1), the first structural 
change at A = Ac witnesses the central vortex being dis- 
placed so as to join the ring formed by the outer ones. A 
further reduction of A leads to a second threshold value, 
\z, at which the mirror symmetry is broken and a zig- 
zag pattern is formed. The situation is different for an 
odd number of vortices, where a parity effect leads to the 
Abrikosov-to-ring and ring-to-zig-zag transitions becom- 
ing degenerate: from full isotropy the lattice re-arranges 
directly into a zig-zag pattern at A = A^, see Fig. [3j For 
even, as well as odd number of vortices a further reduc- 
tion in A makes the vortices align along the weak trapping 
direction, as already observed in Fig. [TJ The situation is 
even richer for a larger (but finite) number of vortices. 
Let us consider, for instance, a system consisting of 18 
vortices (see Fig. |4j). As shown in panel (a), at A = 1 
they arrange in a pattern with a single vortex at the 
centre of the trap and two concentric rings surrounding 
it. By decreasing A we first observe an Abrikosov-to-ring 
structural change involving the inner ring (made out of 
six vortices) and the central vortex, similar to the one 
described above (see panel (b)). By further decreasing 
A, the vortices in the newly formed inner ring start join- 
ing the outer one (see panel (c)) before forming a zig-zag 
pattern (panel (d)). Finally, the transition into a linear 
structure occurs (not shown). 

Let us briefly compare the vortex patterns we have 
just discussed with the ones presented in Refs. [ll|, [l3j . 
In these works the authors address the case of an asym- 
metric trapping potential in the fast rotation limit, Q —> 
mm(uj x , uj y ). In this case the condensate background 
cloud is stretched along the direction of weak confinement 
and assumes a stripe-like shape. Because of the symme- 
try of such a system vortices have to enter the cloud in 
rows and the possible geometries are given by whether 
or not vortices between different rows are aligned. In 
the case of fast rotating traps the solution is determined 
by single particle states and if the vorticity exceeds the 
number of atoms in the system, the existing lattice melts 
and a highly correlated state emerges. In contrast, we 
are dealing with a fixed number of vortices in the limit 



i 
o 




-1 1 
A = 0.70 (c) 




-1 1 
A = 0.3 (d) 




-1 1 



-2 2 



FIG. 4. (Color online) Phase distribution of the condensate 
carrying a vortex- lattice in the x—y plane with N V =1S for 
different values of A. The black dots mark the positions of 
the vortices (in units of y/2N gQ/huo x )- 



where the interaction energy dominates the centrifugal 
one. For such systems the healing length is much smaller 
than any other characteristic length of the system, which, 
as we have shown, leads to a number of possible pat- 
terns with well localised singularities. Transitions be- 
tween these patterns are then determined by the inter- 
play between the trapping potential and the interaction 
energy between the vortices. 



II. SUPERFLUID HYDRODYNAMICS 

In this Section we will explore the structural transi- 
tions in detail by looking at the change in the superfluid 
motion of the condensate. This is analogous to an argu- 
ment used by Fetter in Ref. [HI, where superfluid mo- 
tion in an elliptical and rectangular cylinder was studied. 
While both the e nerg y and angular momentum of the sys- 
tem were found |15j , the existence of a threshold value 
for the angular velocity above which the configuration 
with one vortex is energetically favorable was shown. It 
is important to stress that in our case the vortex-lattice 
configuration found by minimizing Eq. (|2]) does not repre- 
sent, in general, a rigid pattern, due to the perturbations 
introduced into the system by the eccentricity. This can 
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FIG. 5. (Color online) Superfluid velocity field in the rotating 
frame for N v = 8 (other parameters as in Fig.[T]). From (a) to 
(d) the asymmetry parameter is given by A=l, 0.76, 0.56, 0.36. 
Dark purple regions correspond to zero velocity and the ve- 
locities close to the vortex cores are not shown on the chosen 
color- map. 



be seen by recasting the trapping potential as 

V x ( X )^V s ^)+V Q (X,y)=^mu J 2 x (x 2 +y 2 )+^mu J 2 x (X 2 -l)y 2 

(3) 

and recognizing Vq(A, y) as a term exciting quadrupole 
modes. Thus the background condensate and the vortex 
pattern are not stationary. 

The free energy of the rotating BEC is now given 
by Fn v =En v (Q, \)+Ut + Ui where Ut,i are defined by 
Eq. (|2]) and En v (Q.,\) is an energy term that does not 
depend on the vortex configuration and whose detailed 
form is not essential for our discussions. By calling 
{r^} (i = l,..,iVy) the vortex positions which minimize 
Eq. (|2]) for a set number of vortices, we have the condi- 
tion VjF/vJ{ r o\=0, where Vj=(d Xj ,d yj ) and where we 
have used the subscript j to represent the coordinates 
of the j th vortex. In the rotating frame, a vortex has a 
velocity Vj such that 



(4) 



which implies the absence of dissipation, as expected 
from particles moving in a superfluid. A solution 
to this equation is given by v rj =a(VfFN v \ rj ) with 
\/f = (d yj , —d Xj ), where a is the amplitude of the velocity 
field. Its value 
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is found by comparing it with the velocity field 
(h/m)VS— CtXTj in the rotating frame. In this expres- 
sion, S=So+ ^Zi^j @i is the phase of the order parameter 
as seen by the j th vortex, ta,n6j = (y—yj)/(x—Xj) speci- 
fies the polar angle of a reference frame centered on the 
jth vor t ex core [19] an( i 
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is the vortex- free phase of the BEC at position (x,y). 

In Fig. [5] we show the magnitude of the velocity field 
for N v = 8 in a frame which rotates rigidly with the trap. 
The value of A decreases from panel (a) to (d) and the 
arrows show the flow directions with the magnitude being 
encoded in the color. In the dark (dark purple) regions 
the velocity field vanishes, i.e. the superfluid moves at 
the trap angular velocity. For no eccentricity [panel (a)] 
the vortex pattern rotates rigidly with the trap poten- 
tial since the velocity field at the vortex positions (when 
the vortex itself is not present) vanishes in the rotating 
frame. It is worth noticing that outside the vortex pat- 
tern particles flow with a different velocity. This is at the 
origin of the imperfect rigid-body rotation of finite-sized 
vortex patterns in isotropic traps. By increasing the ec- 
centricity [panel (b)-(d)] the rigid body behavior is lost 
and the vortex pattern is no longer a steady solution [tJ , 
since the continuous rotation of the trap increases the 
angular momentum of the system. However, the con- 
dition Q£[Qn v , Qn v +i[ on the angular velocity fixes the 
number of vortices in the condensate N v . The only pos- 
sibility for the system to react is to move the vortex cores 
to accommodate the angular momentum. In a real sys- 
tem, heating and dissipation would eventually lead to the 
crystallization of the vortex pattern or the transition to 
a turbulent regime [22| . 



III. VORTEX LATTICE MODES 

A quantitative confirmation of the abrupt nature of the 
structural changes can be found by studying the eigen- 
modes of the vortex pattern [18| . We take a set of small 
displacements {Sr^ from the equilibrium configuration 
{r^} and write 



6v= (5v x ri , (to* , . . , 5v x , Sv* ) , 



so that the vortex cores velocities in the rotating frame 
become Sv c± A • Sr. Here A is a 2N v x2N v matrix 
whose j th row is found by expanding the velocity field 



N v \t. 



around each r^. This gives 



A^aYJfcUV/FjvJ^ + d Vi (VfF Nv )yi] w , (8) 

i 

where a is determined as before. We now numerically 
diagonalize A for a set number of vortices. The eigen- 
values ai (1<1<2N V ) of A represent the rate at which 
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FIG. 6. (Color online). Spectrum of a BEC with N v = 7 
vortices against the eccentricity A. The points Xc,l where 
the vortex pattern undergoes a structural change are visible. 



vortices start moving from {r^} once they are displaced 
by the corresponding eigenvector 5r l . We note that the 
eigenmodes are related by 

a n (X)+a 2Nv - n (\)=C(X) (0<n<N v ) (9) 

and the corresponding eigenvectors are mutually orthog- 
onal. The constant C(X) depends on the system pa- 
rameters but, remarkably, is independent of the pair of 
eigenvectors considered. A typical spectrum for N v =7 is 
shown in Fig. [6] At two specific values of A the eigen- 
modes show non-continuous behavior, beside the appear- 
ance of a null eigenvalue. These points can be con- 
nected to the structural transition points: Ac signaling 
the Abrikosov-to-ring transition and \l the zig-zag-to- 
linear one. At any other value of A the eigenmodes are 
positive confirming our previous point on the non-steady 
nature of the vortex patterns in the rotating frame. How- 
ever, the exact value of A at which the lowest eigenvalue 
first deviates from zero is found to grow with the number 
of vortices. The corresponding eigenvector corresponds 
to displacements of the vortex positions along the tan- 
gent to the vortex ring, i.e. a rotation of the vortex 
pattern produces no effect. In fact, it is not possible to 
clearly discriminate the eigenmodes of a finite-size lattice 
with a small number of vortices from the phonon modes 
of the background condensate: rotating an anisotropic 
trap excites Bogoliubov modes in the BEC, which have 



a strong influence on the vortex pattern [23|. The link 
between Bogoliubov modes and changes in the proper- 
ties of the vortex matter has already been explored in 
relation to vortex-pattern formation and instability [24j . 
Moreover, in Ref. [25] stability of vortex clusters (com- 
prising both vortices and anti-vortices) in a non rotat- 
ing anisotropic trap has been studied by looking at the 
Bogoliubov modes. The dynamics induced between the 
background cloud and the vortex matter has there been 
shown to be not separable. 

IV. CONCLUSIONS 

We have studied the structural transitions induced in 
a finite vortex-lattice by an increasing degree of eccen- 
tricity of a rotating BEC. An Abrikosov-like arrange- 
ment undergoes a sequence of symmetry-breaking pro- 
cesses that push it towards a linear arrangement of vor- 
tices. Such modifications, witnessed and understood in 
terms of background superfluid motion, are well signaled 
by the eigenmodes of the vortex-lattice. By address- 
ing the case of a finite lattice, our work complements 
and extends the existing literature on vortex instabilities 
and arrangements in rotating BECs and provides inter- 
esting insight into the many-body properties of a meso- 
scopic quantum system. Our analysis is not limited to 
BECs: vortex-like excitations exist in superconducting 
films, Josephson-j unction arrays and dislocation pairs in 
the theory of 2D melting [26|. Inter- vortex potentials 
depending logarithmically on the distance between two 
vortices, similar to Eq. 0, have been observed in thin 
superconducting films [27']. Vortex lattices in thin films 
under magnetic fields have been shown to take the form 
of discrete rows [28|. Strong analogies between the dy- 
namics of vortex lattices and Josephson-j unction arrays 
hold due to the charge- vortex duality [29], thus giving 
our results a generality and interest that goes beyond 
the cases addressed here. 
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